function createFigures()

% 2D pseudo-polar grid for n = 4
% draw 2D pp grid x
name='fig2dx';
figure('name',name);
drawPPgrid(1,0,1,0,0)
save_pic(name);
% draw 2D pp grid y
name='fig2dy';
figure('name',name);
drawPPgrid(1,0,0,1,0)
save_pic(name);
% draw 2D pp grid x and y
name='fig2dxy';
figure('name',name);
drawPPgrid(1,0,1,1,0)
save_pic(name);

% 3D pseudo-polar grid for n = 4
% draw 3D pp grid x
name='fig3dx';
figure('name',name);
drawPPgrid(0,1,1,0,0)
save_pic(name);
% draw 3D pp grid y
name='fig3dy';
figure('name',name);
drawPPgrid(0,1,0,1,0)
save_pic(name);
% draw 3D pp grid z
name='fig3dz';
figure('name',name);
drawPPgrid(0,1,0,0,1)
save_pic(name);
% draw 3D pp grid x, y and z
name='fig3dxyz';
figure('name',name);
drawPPgrid(0,1,1,1,1)
save_pic(name);

%speedup and error with and without preconditioner for 2D
testMakeHFilter2D(1);

%speedup and error with and without preconditioner for 3D
testMakeHFilter3D(1);

%A single cell of the 2D pp grid
drawTrapezoid(2,0);
%A single cell of the 3D pp grid
drawTrapezoid(3,0);

%Calculating the area of a cell in the 2D pp grid
name='calc_trapezoid_area';
figure('name',name);
[I,map] = imread(['latex' filesep 'thesis' filesep name '.png'],'png');
imshow(I,map);
save_pic(name);

%Cells with equal sizes share the same pattern
%Cells of the 2D pp grid
%drawTrapezoid(2,1);
%Cells of the 3D pp grid
%drawTrapezoid(3,1);

%draw all preconditioner related figures
drawPreconditionerRelatedFigures

end
